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Abstract. We calculate the one- loop contributions to the effective potential for the minimal 
■ supersymmetric model when scalar fields other than the Higgses have vacuum expectation values. 

The importance of these contributions for studies of charge and colour breaking bounds is 
discussed. 
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In supersymmetric models, in order to give mass both to the up and down-type quarks (and 
also to cancel chiral and SU(2) anomalies), two Higgs doublets, H\ and H2, are required. The 
neutral components of these fields have a non-zero vev which, much in the same manner as in 
the SM, breaks the electroweak gauge symmetry. However, there are now other scalar fields in 
r S^ • the game, several of them carrying charge, colour or leptonic number, and a priori there is no 

reason why one or many of them cannot have a vev. This problem was first addressed in ref. ffl 
and used to impose restrictions on the undesirably large parameter space of supersymmetric 
models, by demanding that the minimum of the potential along charge and/or colour breaking 
(CCB) directions (meaning, when fields other than H® and acquire vevs) be less deep than 
the so-called "real" minimum. This very appealing idea was developed by many authors 
and of particular importance was the contribution from Gamberini et al ||]: they showed that 
the one-loop corrections to the effective potential in supersymmetric theories were of immense 
relevance for gauge symmetry breaking, in that only at a given renormalisation scale M - namely, 
M of the order of the most significant mass appearing in the one-loop corrections - can their 
contribution to the minimisation conditions be neglected and results derived from the tree-level 
potential trusted. Given that every work on CCB is based on an analysis of the tree-level 
potential only, they concluded that any bound on supersymmetric parameters we obtain in this 
way, if not evaluated at the scale M, will lead to an excessive constriction on the parameter 
space M. A very comprehensive review of the subject is given in ref. M. Recent use of CCB 



bounds in supersymmetric phenomenology may be found in refs. jy, 15]. In this article we will 
investigate the importance of the one-loop contributions to the effective potential in determining 
CCB bounds, as well as discuss the choice of renormalisation scale at which such bounds should 
be considered. 

We will follow the analysis of reference j/J, and use most of the conventions of that article 
(as in the factors of y/2 for the vevs considered). Our convention for the sign of /x is, however, 
the opposite. The superpotential for the model we are considering (no Yukawa couplings other 
than those of the third generation) is given by 

W = X t H 2 Qt R + XbHiQbR + \ t H x Lt r + 11H2H1 

= X t (H+b L - H%t L )t R + \ b (H?b L - H{t L )b R + K(H?t l - H^u l )t r + 



^{H+Hi - Hi Hi) , (1) 
The tree-level effective potential is the sum of the F, D and soft terms, which are given by: 
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Other than the gauginos A a , all fields in the soft potential are scalars. The contributions of the 
first and second generation sparticles are not zero - for simplicity we do not write them, but they 
are identical in form to those of their third generation counterparts, minus all terms proportional 
to Yukawa couplings, which are set to zero. The one-loop contributions, at a renormalisation 
scale M, are given by 



where the M a are the (tree-level) masses of each particle of spin s a and n a = (— l) 2sa (2s a + 
1) C a Q a . C a is the number of colour degrees of freedom (1 for colourless particles, 3 for all 
others except the gluinos, which, being an SU(3)c octet, correspond to a factor of 8) and Q a is 
2 for charged particles, 1 for chargeless ones. Thus the factor n a accounts for the multiplicity 
of each particle. When the neutral components of H\ and H 2 acquire vevs vi/y/2 and v 2 /\< r 2, 
we obtain the "real" minimum, and the tree-level potential Vq may be written as 



t r 1 / 2 2, 22 

= 2 V 1 Vl m2 V2 



Bfiviv 2 + —{g 2 + gl) {vl - vl) 2 



(6) 



with m 2 = m 2 Hi + fj? and m 2 = m 2 H2 + fj?, and where the vevs' phases were chosen such that the 
term with the bilinear coupling is negative. Minimising with respect to v\ and v 2 we obtain, for 



the value of the potential at the minimum, 



v^ SM = -JLtf 2 + pi) («£-«?)* . (7) 

So that a CCB minimum can exist, we need negative terms in the potential other than the one 
containing B \i. These arise from trilinear terms in the potential and require non-zero vevs for 
at least two fields other than H® and H®. As a working example, consider the case where the 
fields tl and tr have vevs (l/y/2 and r/y/2 respectively). The vacuum tree- level potential is 
now given by 

Vo = %vI(1 2 + t 2 )+1 2 t 2 } - ^(A TVl + tiv 2 )lT + \{m\v\ + m 2 v 2 + m\l 2 + 

m 2 T r 2 ) - B„ Vl v 2 + 9 ^(v 2 -v 2 -l 2 + 2r 2 ) 2 + | (v 2 - v 2 + I 2 ) 2 , (8) 

where we chose the phases of the vevs so that the terms with bilinear and trilinear couplings 
are negative []. With four different vevs in the potential, an analytical study of CCB is very 
difficult if not impossible. Following || we express all vevs in terms of, say, v\, so that I = a v\, 
V2 = 0v% and t = jv±, with constant and positive parameters a, j3, 7. The potential ||] becomes 
a quartic polynomial in v\ and its minimisation a trivial exercise. Comparing the value of 
the CCB minimum with that of eq. ^ and requiring this to be smaller, we find bounds on the 
MSSM parameters as function of a, . . . p|. The vevs are found to be of order v ~ |A T |/A T , 
and the potential V ~ — | Aj-^/A 2 .. An important observation throughout is the choice of the 
renormalisation scale at which this calculation is made: as first shown in Q, the tree-level 
vevs are only reliable at a renormalisation scale of the order of the largest mass present in 
the one-loop contributions |5[ For our CCB potential, we expect this to be of order M ~ 
max(At , g2)|A T |/4A T ||. However, the highest mass present in the MSSM potential will be of 
the order of the gaugino mass, Mq, in principle different from M. This line of thought therefore 
dictates the two potentials should be compared at different renormalisation scales. Now, the 
work of Gamberini et al refers to the MSSM and demonstrates quite clearly that the tree-level 
vevs can be trusted in a given range of the renormalisation scale. Meaning, that in that range 
the one- loop contributions to the minimisation conditions, to wit 

Y^M 2 ^(lo g ^-l) (9) 
% 32vr 2 a dvi { g M 2 J ' [ } 

are not significant. It does not, however, say anything about the size of the one-loop contributions 
to the potential themselves. We expect them to be small, given that the choice of M limits the 
size of the logarithmic terms, but notice this is not guaranteed. Another and stronger argument 



for the inclusion of the one-loop potential in CCB calculations comes from ref. [11|. Their results 
are thus summarised: the sum V0+AV1 is not one-loop renormalisation group invariant. Instead, 
the full RGE invariant effective potential is given by 

ViM^fa) = n(M,Xi) + VoiXi,^) + AVi(M,Xi,(j)j) + 0(h 2 ) , (10) 

where Aj stands for all couplings and masses of the theory and (f)j for all its fields. The function 
depends in principle on M - implicitly or explicitly - and Aj, and receives contributions from 
all orders of perturbation theory. It is, however, field-independent 0. The authors of 0] give 

Depending on the signs of A T , fi and B, this isn't always possible. This is however a detail, and it doesn't 
have any impact on the conclusions we will draw. 

2 Which of course means it usually has no effect whatsoever in particle phenomenology, given all quantities of 
physical interest are given by derivatives of V with respect to the fields. 



examples of explicit calculations of f2 for different models and its importance. Notice now that 
the only difference between the "real" MSSM potential and the CCB one is that some fields have 
vevs, others do not - in other words, the two potentials correspond to two different sets of values 
for the fields in the theory. But the function f2 is field-independent and is thus the same, CCB 
or not. So, when we compare y MSSM and V CCB at different renormalisation scales, we are 
not taking into account the contribution arising from Q. We are led to the conclusion that, to 
eliminate f2 altogether from the calculation, we must compare 

yMSSM and yCCB ^ thfi 

same 



renormalisation scale. Also, notice that eq. 10 ensures that if y MSSM > yCCB f or a gi ven sca l e 
M, the inequality holds for all scales. Now, the rationale we gave for comparing both potentials 
at different scales was the possibility of neglecting the complex one- loop contributions. The 
existence of forcing a comparison at equal M, we conclude that, for one of the potentials 
at least, the full one-loop vevs must be used. By consistency we should then work with the 
value of the potential correct to one-loop as well. Unfortunately, this way we abandon all hope 
of obtaining analytical CCB bounds, given the complexity of the one-loop contributions. We 
now present their calculation in the CCB case. Non-zero vevs will be considered only for the 
real components of the fields Hf (v\/V2), H® (V2/V2), tl (1/V2) and tr (t/V2). This choice 
should be favourable to CCB in that and m T are generally the smallest of the soft masses 
and, because no squarks have vevs, there are no SU(3) D-terms, thus limiting the positive 
contributions to the potential. Also, A T being typically quite small, a CCB minimum, should 
it exist, ought to be quite deep. The simplest sparticle mass is that of the gluinos (n g = —16), 
it is simply M3, from eq. ||. The gluino contribution to V\ is thus field-independent and could 
be included in the fi function in eq. |l0| - but in any case where colour breaking does occur 
the gluinos would mix with other fields, so, we prefer to single out this contribution 3 . Also 
very simple are the contributions from the top and bottom quarks (rit = 11^ = —12), with their 
squared masses given, as usual, by A 2 v 2 /2 and \\v\j2 respectively. For the stop (n-i = n% = 6), 
we have 
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The equality £12 = £21 arises from the fact we are considering all parameters in the potential to 
be real (the same happens for all the other mass matrices, who are thus symmetric). For the 
sbottom (rii = 112 = 6), 
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3 Of course, given that y CCB and y MSSM w i\[ be compared at the same renormalisation scale, the gluino 
contribution will be identical in both potentials and have no impact on CCB. 
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The masses for the first and second generations squarks are obtained from these matrices by 
setting the Yukawas to zero and replacing the soft masses by the appropriate ones. The squarks 
are degenerate for these two generations (assuming universality of the soft masses - but these 
results can be trivially generalised for non-universal models) and, for the up-type squarks (ni = 
n2 = 2 x 6 - the extra factor of two accounts for both generations), we have 
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Likewise the first and second generation sleptons (ni = n 2 = 2 x 2) will have degenerate masses 
given by 



Ml 
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In these expressions, mR, m u , m^, mjy and m e are the lower-generation equivalents of ttiq, mj, 
m^, and m T respectively. Finally, the first and second generation sneutrinos (n = 2 x 2), 
also degenerate, have a mass given by 



Mr, 
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The masses of the gauge bosons are easily determined (see, for example, |T^]). For the charged 
gauge bosons (we still call them W's; n = 6) we have 



Mi 
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The mass matrix of the neutral gauge bosons (n\ = n 2 = 3) is given by 
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Notice that if I = r = this matrix produces the usual photon and Z masses. The remaining 
sparticles' mass matrices are complicated by mixing between fields. Even when there is no CCB 
some mixing occurs - the charginos are after all the result of the mixing between the charged 
gauginos and fermionic partners of the Higgses, the neutralinos a similar mix but between neutral 
fields - but with the vevs I and r in the game, mixing between neutral and charged particles 
occurs. In the case of the charginos, for instance, the interaction term in the Lagrangian between 
the gauge and matter multiplets, plus the Yukawa interactions in the superpotential, cause a 



mixing with the r neutrino (see, for example, [13]), vl (we use the bar so as to not confuse this 
field with ul, one of the components of the scalar doublet L) so that, if X^ are the charged 
gauginos and the fermionic partners of the charged Higgs fields, the chargino mass matrix 
will be given by (the columns, from left to right, correspond to the fields A H 
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This matrix has one zero eigenvalue (corresponding to the neutrino which, despite the mixing, 
does not gain a mass in this scenario) and the squared masses of the two charginos (rti = n>2 = 
—4) are determined by a simple quadratic equation. The neutralinos (ni...6 = —2) now include 
mixing with the tau fermion, and their mass matrix is (from left to right, the entries of the 
matrix are B, W 3 , Hi 
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The charged Higgs (ni...3 = 2) fields of the theory mix with the r sneutrinos, so that their 
combined mass matrix is (from left to right, H r , H 2 and ul) 
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It is easy to check, using the tree-level minimisation conditions, that this matrix has a zero 
eigenvalue, corresponding to the two Goldstone bosons that give mass to the charged W gauge 
bosons. The pseudoscalars (ni...4 = 1) in this theory are a mix between the imaginary parts of 
the neutral components of Hi, Hi and tl and tr. The mass matrix is given by 
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where each of the entries of this matrix is a derivative of the form d 2 V/dIm i 
fields (pi corresponding to, from left to right, H®, H®, tl and tr. We have 
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Using the tree- level minimisation conditions, we see that this matrix has two zero eigenvalues, 
corresponding to the Goldstone bosons that give mass to the neutral gauge bosons - we thus 
have four Goldstone bosons, one charged and two neutral, exactly as we should have, given we 
are fully breaking the SU(2)\y x U(l)y gauge symmetry. The "Higgs scalars" (ni„A = 1) of 
this theory are a mix of the real parts of the neutral components of Hi and H 2 and tl and tr. 
The mass matrix is 
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dn Ih h H ih 
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A simple check we can perform on the validity of these formulas is to calculate Str M 2 = 
J2 a n a M 2 - because supersymmetry is broken in a "soft" way, this quantity should be field- 
independent. It is not difficult to verify that the mass matrices here presented satisfy this 
condition. Let us now apply these results to a simple example, the MSSM with universality 
(common values for the gaugino and scalar masses and the A parameters at the gauge unifica- 
tion scale My ~ 10 16 GeV). Besides their simplicity there are several reasons these models are of 
interest: the absence of flavour changing neutral currents at the weak scale is a strong indicator 
of universal soft parameters. There is also the evidence that a vast range of theories have an 
universal form at My as an infrared fixed point []l4|]. Finally, these models are of phenomeno- 
logical interest in that their predictions for the sparticle masses have not been excluded by the 
current experimental data (see, for instance, ref. []l^]). A brief description of our procedure: 
we input at scale Mz the masses of the gauge bosons, third generation fermions, the values of 
the gauge couplings and tan/3 (ranging from 2.5 to 10.5). Using the two-loop /3-functions, we 
then find the gauge unification scale (unification of a\ and 02 only) at which point we input 
the universal soft parameters (20 < M G < 100 GeV, 10 < m G < 160 GeV, -600 < A G < 600 
GeV). We then go down to scale Ms = max(M^, M G , rn G ) (of the order of the largest mass in 
the theory) and determine B and fi from the minimisation of the full one-loop potential (see, 
for instance, (jTj). At this stage we have about 3200 "points" of the parameter space of the 
MSSM. At the scale M = 0.6 g2\A T \/X T we then solve the CCB minimisation conditions - the 
set of equations {dVo/dvi = , dVo/dv2 = , dVo/dl = , dVo/dr = 0}, with Vq from eq. || 
- which is equivalent to determining the vevs {t>i , v 2 , / , r}. Notice that we solve the tree-level 
minimisation conditions for the CCB potential - we are relying on the work of Gamberini et 
al, assuming that at the scale M the contributions |9| are not significative 0. We then calculate 
(Vq + AVi) CCB and compare it - at the same scale M - with (Vq + AV\) MSSM . For much of the 
parameter space considered solving the CCB minimisation conditions is not possible. For ~ 1290 
"points", we find extrema for which the value of the potential is smaller than (Vq + AV\) M . 
We can evaluate the importance the choice of renormalisation scale has - if we choose, for in- 
stance, M = 52 1 A T I /4A T , we see CCB occurring for 200 more "points". Notice how our choice 
of renormalisation scale is always a field-independent one This dependence on M would, 

4 We chose M, after several attempts, to be the value for which the ratio r — max(Afa)/M 2 is closest to 1. At 
this scale the one-loop contributions to the minimisation of y MSSM are expected to be significative, so, they are 



included. 



o 



one expects, vanish if we had performed a full one-loop minimisation of the CCB potential [|]. 

These CCB extrema are distributed uniformly through the input values of tan (5, Mq and mc, 
but concentrate on the high absolute values of Aq - we only see then occurring for \Aq\ > 300 
GeV. This of course reminds us of the "classic" CCB bound according to which CCB would 
occur if \A T \ 2 > 3 (m\ + m\ + m^). Figure p] shows us that the majority of rejected "points" 
doesn't satisfy this bound. They also concentrate in specific values of the \x parameter, as may 
be seen in fig. |2[ Interestingly, the rejected "points" target higher values of - notice how 
few rejections occur for Mh < 80 GeV - and do not occur for small values of \fi\. This Higgs 
mass there plotted is that of the "real" , MSSM, Higgs boson, calculated from the full one-loop 
mass matrices. Another interesting observation is that the values found for the vevs I and r are 
very similar. To evaluate the impact of using the full one-loop potential and comparing both 
potentials at the same scale we performed the "usual" CCB analysis (i.e, compare the tree-level 
potential only, at two different renormalisation scales) as well. We found CCB happening for ~ 
10% more "points" than we have found using the full one-loop potential. 

Some remarks: the CCB extrema we found are not minima, in that there are always some 
quadratic sparticle masses which are negative (specifically, &2> the m's and v e ). It has been 
shown that a non-convex effective potential is the result of a bad summation of 1PI diagrams, 
and the way to proceed in such situations is to take the convex envelope of the potential. A 
good approximation is to consider the potential flat where negative second derivatives occur. We 
point out, however, that the negative masses found are very small when compared to the typical 
masses in the CCB potential, and can therefore be safely neglected. Also, we must remember 
that all that has been said here relates to the possibility of the "real" vacuum being the absolute 
minimum of the potential. There is however the argument |2| that CCB bounds should be 
determined from requiring the tunnelling time from the "real" minimum to the CCB one to be 
larger than the age of the universe. This, of course, leads to a relaxation on the bounds. We 
do not undertake this study here, but argue that, in any case, the one-loop contributions to the 
potential would be relevant, in that they determine how deep the CCB minima/extrema are. The 
importance of the one-loop contributions to CCB bounds had already been studied in ref. [J16[|, 
but comparison with their results is difficult: they considered only the top-stop contributions 
to AVi (we found that those are very important, but no more so than the chargino and "Higgs 
scalars" ones) and used a calculational method very different from the one here presented. In 
conclusion, we calculated a full one-loop CCB potential, and argued that the existence of the 
function - that is to say, the RG invariance of the effective potential - forces a comparison of 
the CCB and MSSM potentials at the same renormalisation scale M. Having the actual CCB 
masses informed our choice of M, higher than usually considered. In a simple example, we found 
a small, but nevertheless sizeable, difference between the parameter space rejected in the two 
CCB approaches here delineated. It remains to be seen what that difference might be for other 
models, more complex or with larger parameter space. In light of the (hopefully) impending 
discovery of the Higgs boson, a careful reconsideration of CCB bounds might be in order. Finally, 
the computation of the CCB mass matrices opens the way to the possibility of performing a full 
one-loop minimisation of the CCB potential. That way one could confirm if their impact in the 
bounds is indeed negligible, and eliminate the renormalisation scale dependence we encountered. 
This work is in progress 
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Figure Captions. 

Fig. [l]. A 2 r versus Q = 3(m± + m 2 L + m 2 .), for the rejected CCB points at the scale M; all 
points below the solid line violate the "classic" CCB bound of Frere et al. 

Fig. [2]. Mass of the lightest Higgs boson versus the value of the \i parameter at the unification 
scale, for the CCB-rejected points. 
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Figure 1: A 2 T versus Q = 3 {m\ + m? L + m 2 ), for the rejected CCB points at the scale M; all 
points below the solid line violate the "classic" CCB bound of Frere et al. 




Figure 2: Mass of the lightest Higgs boson versus the value of the [i parameter at the unification 
scale, for the CCB-rejected points. 
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